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I. INTRODUCTION 

Since the pioneering papers of Polyakov [1] and Susskind 
[2] the confinement-deconfinement phase transition in finite 
temperature Yang-Mills theory has become a thoroughly stud- 
ied phenomenon. This is particularly true for the 'standard' 
groups SU (2) and SU (3), but more 'exotic' groups like G(2) 
have recently come into focus as well, see e.g. [3, 4]. Due to 
the nonperturbative nature of the problem progress has mainly 
been achieved via brute force lattice computations. 

Nevertheless, it would be helpful to have a simpler and 
more intuitive understanding of the physics involved. The 
principal tool for this purpose is the construction and subse- 
quent analysis of effective models. Such attempts also have 
quite some history going back to the works of Svetitsky and 
Yaffe [5, 6] as well as Polonyi and Szlachanyi [7]. The basic 
idea (in the spirit of Landau and Ginzburg) is to use the order 
parameter of the transition, the Polyakov loop, as a collective 
degree of freedom and formulate effective theories in terms of 
it. For gauge groups SU (N) the rationale behind that is the 
Svetitsky-Yaffe conjecture [5, 6] which states that the Yang- 
Mills finite-temperature transition in dimension d + 1 is de- 
scribed by an effective spin model in d dimensions with short- 
range interactions 1 . These ideas have initially been taken up in 
terms of strong coupling expansions [7, 9, 10] yielding Ising 
type spin models with an effective coupling A(/3) where (3 
denotes the Yang-Mills-Wilson coupling. A review of early 
work in this context may be found in [1 1]. For an overview of 
more recent developments we refer the reader to [12]. 

Early on, it has also been attempted to obtain these ef- 
fective models, that is their couplings (being the 'weights' 
of the included operators) nonperturbatively via lattice meth- 
ods. In [13, 14] Creutz's microcanonical demon method [15] 
has been employed for SU(2). An alternative method based 



The reasoning involved strongly relies on center symmetry. The study of 
more exotic Lie groups (which may not even have a nontrivial center) sug- 
gests that the size of the gauge group is also important [8] . 



on Schwinger-Dyson equations (SDEs) and dubbed "inverse 
Monte Carlo" (IMC) was developed soon after [16, 17] and 
applied to both SU{2) [18, 19] and 5/7(3) [20, 21]. Since 
then the IMC approach to lattice gauge theories has largely 
been dormant with only a few exceptions [22, 23]. 

Inspired by the success of Polyakov loop models [24-26] 
we have recently reinvestigated the feasibility of IMC for the 
confinement-deconfinement phase transition in a series of pa- 
pers. Our numerical approach has consistently been comple- 
mented by analytical attempts such as strong coupling expan- 
sions and mean-field approximations. For the second-order 
SU(2) transition our results may be found in [27] and [28]. 
By including up to 14 operators and 3 different group rep- 
resentations we were able to reproduce suitable Yang-Mills 
observables to a reasonable accuracy. The same is true for 
the analytically known asymptotic behavior of the effective 
couplings as a function of (3. In [29] we have started to inves- 
tigate effective models for SU(3) which are generalisations of 
the 3-states Potts model. The critical behavior of these is an 
interesting subject in its own right. We have found a very rich 
phase structure with first and second order transitions between 
symmetric, ferromagnetic and anti-ferromagnetic phases. In 
addition there seems to be a tricritical point rendering mean- 
field theory approximately exact in its vicinity [29, 30] 

In relating the effective models to SU (3) Yang-Mills via 
IMC one expects to encounter new difficulties. The first tech- 
nical problem to overcome is to find the SDEs which are less 
straightforward than for SU (2) as the SU(3) group manifold 
no longer is a sphere. This problem has been solved in [31] 
and [32]. As the SU(3) phase transition is of (weak) first or- 
der, hence not continuous, the determination of the effective 
couplings might require fine-tuning raising the question of sta- 
bility of the solutions. This will be one of the main issues to 
be addressed in what follows. 

The remainder of the paper is organized as follows. In Sec- 
tion II we suggest different effective actions as candidates for 
describing the Polyakov loop dynamics of Yang-Mills theory. 
In Section III we explain how to obtain the effective couplings 
via two alternative sets of SDEs and subsequent IMC. Our nu- 
merical results are presented in Section IV. Section V con- 
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eludes our discussion with a summary and outlook. Some 
technicalities are deferred to Appendices A-C. 



II. EFFECTIVE ACTIONS 

We begin by recalling the lattice definition of the untraced 
Polyakov loop in the group representation 7Z, 

nw x ) = Y[K(u tx ) , (i) 

where 7Z(Ut x ) is the temporal link at time slice t and posi- 
tion x in representation 7Z. Any irreducible representation 
of 5/7(3) is labeled by two integers, 1Z = lZ pq , which, in 
flavor language, count the number of quarks and antiquarks 
needed to construct the multiplet associated with TZ pq . The 
basic building blocks of our effective actions are the group 
characters xn associated with the representation 1Z, that is 
the traces of the Polyakov loop (1), 

X n(V) = x Pq (V)=trK pq (¥) . (2) 

Note that these only depend on the traced Polyakov loop in the 
fundamental representation, V = tr *p. The trivial character is 
of course xoo = 1 while the first nontrivial ones correspond to 
the (anti-)fundamental representations and yield just the stan- 
dard traced Polyakov loop (and its complex conjugate), 

Xio(V)=V, Xoi(V)=V*. (3) 
Under a center transformation, the characters transform as 

Xpq Z^Xpq , ZtZ 3 . (4) 

Center symmetry is then sufficient to determine the opera- 
tor content of the effective action if we restrict to nearest- 
neighbor (NN) interactions. In terms of group characters one 
finds the general form 

,p'q' Xpq {V^Xp'q'iVy) , (5) 

(xy),pq,p'q' 
p+p' —q-\-q' mod 3 

where the sum over representations is constrained by center 
symmetry. Expressing the characters explicitly in terms of 
the Polyakov loop V one easily recognizes (5) as the action 
suggested by Dumitru et al. [24]. Their 'potential terms', 
built from single center symmetric characters located at single 
sites appear whenever the second adjacent character is trivial, 
Xoo = 1. In this case the typical hopping terms connecting 
NN sites 'degenerate' into ultra-local terms, the one of lowest 
dimension being the 'octet loop' contribution, Au^o Xn as 
indeed 1 + = 1 + mod 3. One expects that the couplings 
\ P q,p>q> decrease with increasing representation labels, p, q, p' 
and q' , hence that representations of low dimension, 

d pq = \{p+l)(q+l){p + q + 2) , (6) 



dominate the effective action. To simplify our notation we will 
henceforth write the action (5) (and generalizations thereof) as 
a series of the form 

S e ff = X Xi Si , (7) 

i 

where up to 16 different terms Si will be taken into consider- 
ation, albeit not necessarily within one and the same ansatz. 
A list of the operators Si may be found in Appendix A where 
we allow for next-to-NN (NNN) couplings in addition. It is 
easy to check that each of the terms given satisfies the se- 
lection rules for the representation labels necessary for center 
symmetry. 

It turns out (in hindsight) that the ansatze (5) or (7) contain 
more freedom of choice than actually required which makes 
the inverse Monte Carlo routines less efficient (see below). 
To further constrain this freedom we generalise to SU(3) 
the strong-coupling approach introduced by Billo et al. [33] 
which we already have successfully applied to SU(2) [28]. 
The basic building blocks are then given by the center sym- 
metric operators connecting NN sites, 

S n ,i = XTz(V x )x*Tz(V y ) + c.c. , £ = (xy) . (8) 

The strong coupling expansion then replaces the ansatz (5) by 
the following somewhat more complicated expression [31], 

r 

£* = £ £ £ 4;:^ r (/3)II^A, (9) 

r K-L...K r li—Zr i=l 

where r counts the number of link operators (8) contribut- 
ing at each order. The coefficients c^' " lr n are the couplings 
between the operators S-n i ^ i from (8) sitting at NN links 
li = (xi,yi) in representation IZi. The effective action de- 
fined in (9) hence describes a network of link operators of 
the type (8) that are collected into (possibly disconnected) 
'polymers' contributing with 'weight' ^ . Again, one 
expects the 'weights' or couplings to decrease as the dimen- 
sions of the representations and the inter-link distances in- 
volved increase. In a strong coupling (small-/?) expansion 
truncated at 0(j3 kNt ) one has r < k and the additional re- 
striction \1Zi | + --- + \lZ r | < k with \1Z\ = p + q for a given 
representation 7?.. 

To lowest order [3 Nt one finds the universal effective action 
[7] 

S e n = Cio X S W ,(xy) = «1 ^2(Xlo(Px)XOlCPy)+£-£-) 
(xy) (xy) 

(xy) 

which is just a single hopping term connecting Polyakov loops 
at NN sites. This is reminiscent of a generalised Ising model 
or, more appropriately, a three-state Potts model [34]. As 
mentioned before the study of these models is interesting in 
its own right [29] but will not be the topic of the present paper 
which focuses on the relation between the effective actions 
and Yang-Mills theory. 
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Again, the most general representation (9) is not too illumi- 
nating and we will therefore adopt the notation 



5eff 



(ID 



A list of the leading NN and NNN action terms I a can be 
found in Appendix B. 

In principle, without any truncations, the effective actions 
(7) and (11) have to coincide although the operator bases used 
are different. Accordingly, there is a linear relationship be- 
tween the couplings A; and n a . However, as the ordering 
principles for the two ansatze are not the same the relation 
between the couplings after truncation is not one-to-one but 
rather of the form 



A, — K;, 



(12) 



Typically, for a given truncation of the A-action (7) the range 
of the index i is larger than that of a, i.e. there are more A's 
than re's. Accordingly, the matrices (Kia) are rectangular with 
a < i and integer entries. For this reason it is calculationally 
often more efficient to work with the K-action (11) and reob- 
tain the A^ via (12). Our standard choices for the numerical 
matrices (Ki a ) corresponding to different truncations of the 
effective actions may be found in Appendix C. As the opera- 
tors Si appearing in (7) are more intuitive and resemble gen- 
eralised spin terms we have decided, for the sake of brevity, to 
present results only for the A's in this paper. 



III. SCHWINGER-DYSON EQUATIONS AND INVERSE 
MONTE-CARLO 

In this section we shortly recapitulate the SU(3) 
Schwinger-Dyson equations (SDEs) that have recently been 
derived in [31] and [32]. They will be the main tool to re- 
late the effective actions (7) and (11) to Yang-Mills theory. 
Our numerical approach benefits from the fact that there are 
two independent versions of SDEs which in the end, however, 
should yield equivalent results. The first type of equations is 
based on an integral identity which is more algebraic in nature 
than the second type which follows from geometrical consid- 
erations. 



A. Algebraic SDEs 

It is useful to parameterize the diagonalized, untraced 
Polyakov loop by means of two angular variables, (pi and <p2, 



/ e i<t>i o 

Wfafa) = ( 

e-Hfa+to) 



(13) 



with values in a fundamental region given by the restrictions 

4>i < 4>2 < [~4>i - 4>2) mod 2tt , < tpi < 2?r . (14) 



As a result, the reduced Haar measure acquires the form 

dfired = J 2 d(/)id4>2 (15) 

with a nontrivial Jacobian that may either be expressed in 
terms of characters, 



J =15- 6xn + 3x30 + 3xo3 ~ X22 , 
or in terms of the trace V of (13), 

J 2 = 27 — 18VV* + 4"P 3 + 4V* 3 — V 2 V* 2 . 



(16) 



(17) 



Using the latter variable leads to the remarkable algebraic 
identity 

dPredfai.tfa) = J 2 d4>i dfo = JdVdV* = d» ieA {V,V*) , 

(18) 

such that one can trade J 2 for its square root J. It is a fact of 
life that any function / = f(V,V*) vanishing on the bound- 
ary dil of a region SI satisfies the integral identity 



dVdV*d v f = 



(19) 



which is reminiscent of integration by parts on the real line, 
/ dx f'(x) = (for functions / vanishing at infinity). For our 
purposes we make the particular choice 



f(V,V*) = J 3 g{V,V* 



(20) 



with arbitrary g and integrate over the domain of V given im- 
plicitly via (14). This is consistent with the general identity 
(19) as the reduced Haar measure vanishes at the boundary 
<9£! of the fundamental region (14). Hence, (19) specializes to 







J dfi Ied (V z ,V* z 



3 dJl 



2dV 2 



.9 + J; 



dg_ 



(21) 



where we have used the reduced Haar measure from (18) and 
reinstated the dependence on the lattice site chosen to be z. 
The derivative of J 2 can actually be worked out with the result 



dJ 2 
dV 



12V 2 - 2VV* 2 - 18V* 



(22) 



To actually obtain genuine SDEs we have to introduce the 
usual Boltzmann factor exp(— 5 e ff). This is done by choos- 
ing a special function g, namely 



3a 



dh 
dVZ 



exp(-Seff) , 



(23) 



with another function h(V, V*) required to be Z3 invariant, 

h(V,V*) =h(zV,z*V*) , zeZ 3 . (24) 

Plugging in the ansatz (7) for the A-action the P-derivative of 
(23) needed for (21) becomes 



dg x 
dV z 



-Seft 



(25) 
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where the commas denote differentiation with respect to the 
subsequent argument. Functional integration of (21) with the 
measure 3t(i lc d = Y\ z d^ lc d{'Pz,'Pz) finally yields the de- 
sired SDEs, 

/ o a r2 \ 

^ ' i 

(26) 

which comprise a linear system for the effective couplings Xi, 
generalising analogous results for SU(2) [27, 28]. The expe- 
rience gained there prompts us to choose the function h from 
the operators Si in the ansatz for the A-action (7). Any index 
i then yields an independent equation. In addition, this choice 
automatically satisfies the criterion (24) of Z3 invariance. 

On top of that we will vary the sites x and z, in par- 
ticular the distance d = \x — z\ between them. On a lat- 
tice with spatial extent N s this implies a range of distances 
d € {0, . . . , [jV s / 2 J } where [x\ denotes the largest integer 
< x. For TV different operators S{ we thus obtain N indepen- 
dent equations for each distance d. 



which will serve as a master identity generating all SDEs. As 
in the previous subsection we would like to integrate over the 
reduced Haar measure d/i rc d only. Thus, we want the inte- 
grand L a f(U) to be a class function. If G is such a class 
function it only depends on the fundamental group characters, 
Xf(U) = tr(72.i? ([/)), where 1Zf(U) denotes a fundamental 
representation of the group element U, i.e. 3 or 3 for SU(3). 
For the particular choice / = gL a x P , with g an arbitrary class 
function and Xp a fundamental character, the integrand in (28) 
indeed becomes a class function so that we end up with the 
'reduced' integral 

J d^ cd L a {gL a Xp )^0 . (29) 

For SU (3) we obviously choose the fundamental character Xp 
as the trace of the Polyakov loop in the fundamental represen- 
tation, Xp = X10 = ~P an d g = h e~ SdT , slightly different from 
(23). The left derivatives are worked out as follows [31, 32], 



B. Geometrical SDEs 

For any function / = f(U) on a Lie group we define its left 
derivative in the direction of the generator T a via 



L a f(U) = -f (e tJ U 



t=o 



(27) 



Left invariance of the Haar measure implies a symmetry rela- 
tion somewhat analogous to (19), 



d maai {U)L a f(U) = Q , 



(28) 



L a (L a (v)) = -fr , 

L a (V)L a (V) = AV* - ^r 2 , 
L a (V)L a (V*) = 6~l\V\ 2 , 



(30) 



and only depend on the traced Polyakov loop as they should. 
Finally, to obtain feasible equations we choose h among the 
V -derivatives of the Si, h = Si-p implying the following set 
of geometrical SDEs, 



0= - 



Y^ s i,r. + m - ^DS^.v* + (6 - % \7 
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{*K - iPDSi.v.S^v, + (6 - - iVzHSi^Sj^ 



(31) 



where, again, the dependence on the lattice site z has been 
made explicit. 



C. Normalisation 

We have seen that, for every pair of lattice sites x and y, we 
end up with a linear system of equations for the couplings of 
the effective theory. Since on the lattice we have both trans- 
lational and (discrete) rotational symmetry it is sufficient to 
consider different distances d = \x — y\ only. These serve as 
a label for our sets of equations which in a condensed matrix 



notation may be written as 

M d X = b d . (32) 

If we assume a total of N unknown couplings collected into 
the vector A we have, by construction of the SDEs, an N x N 
matrix Md and an inhomogeneity bd, hence an independent 
system of equations, for each distance d. The off-diagonal 
entries of Md and the vector bd are typically complex but the 
couplings A have to be real. We therefore distinguish between 
real and imaginary parts of (32) for different d and group them 



5 



together into the equations 

/ Re M \ 
Im M 
Re Mi 
ImMi 



\Im Mi 



LiV s /2j 



Re6 
Im bo 
Re bi 
Im b\ 

\Imb L jv s /2j/ 



(33) 



This constitutes an overdetermined linear system of 2N x 
{\_N s /2\ + 1) equations for the N unknown couplings A. In 
principle, this can be solved by standard least-square methods. 

However, this procedure is hampered by a few techni- 
cal pitfalls. Since the order parameter for the confinement- 
deconfinement transition is driven by the long-range behavior 
of the lattice system we have to take into account the fact that 
equations associated with different distances d enter (33) with 
different multiplicities. On a three-dimensional lattice this en- 
tails that an equation for distance d > has multiplicity 



m d = (2d+l) 3 - (2d - l) a 



(34) 



while equations for d = only appear once. To account for 
this mismatch we reweight the coefficients of bd and Md (for 
d > 0) with a factor y'mj. 

Another problem are the large condition numbers of the ma- 
trices. To cope with this we employ a simple form of normal- 
ization based on the diagonal elements of the matrix Mo. For 
the algebraic SDEs the diagonal entries Mo,m of Mo dominate 
the linear system. For this reason we construct a new diagonal 
matrix Ao from the Mo.™ according to 

A = diag(Re (M , u )- 1/2 , . . . ,Re (M ,nn)- 1/2 ) ■ (35) 

By means of a similarity transformation with Nq the real parts 
of the Mo,a may be transformed to unity. As a result (33) 
becomes 



( A M d N ) (Nq 1 A) = (N M d N ) ^ = N b d 



(36) 



Typically, this new system of equations is better conditioned 
which increases the stability of the results. 

In our numerical calculations we have both used the im- 
provement (36) in condition numbers and the reweighting fac- 
tors given by the square root of (34). 



IV. NUMERICAL RESULTS 

The IMC method basically amounts to solving the SDEs 
(26) or (31), the crux being the evaluation of the expectation 
values (. . .) in the microscopic ensemble. In our case this is 
given by a sufficient number of SU (3) Yang-Mills configu- 
rations generated by standard MC routines [35-37]. In what 
follows we will consider three ansatze, one with five NN cou- 
plings and two more general ones which either contain NN 
terms in higher representations or additional NNN couplings. 
In a strong coupling expansion (which is rationale for the k- 
actions) both the five leading NN terms as well as the NNN 
ones would be of order (3 2Nt while the extended NN ones 
would beO(P 3Nt ). 



A. NN couplings 

Before one determines the effective couplings correspond- 
ing to Yang-Mills it is prudent to check if the SDEs (26) and 
(31) derived above are consistent within the effective theories 
themselves. To test for that we have first simulated an effec- 
tive theory with five fixed input couplings and tried to repro- 
duce them via IMC. In Table I we have listed the outcome of 
this testing procedure. One notes that the couplings Ajjmc 
determined via IMC coincide with the chosen input couplings 
Ai, input to an accuracy of about 2%, both for the algebraic 
and geometrical SDEs. This tells us two things, first that our 
SDEs (26) and (31) are both correct and, second, that IMC 
works extremely well for the effective theories. 



TABLE I: IMC consistency check for the NN A-action including five 
couplings Ai, . . . , A5. 



algebraic 



geometrical 



i 


Aijinput 


Ai,IMC 


At, input 


A;, IMC 


1 


-0.0100 


-0.0101(1) 


-0.0100 


-0.0101(1) 


2 


0.0060 


0.0059(1) 


0.0060 


0.0060(1) 


3 


-0.0050 


-0.0051(1) 


-0.0050 


-0.0051(1) 


4 


0.0080 


0.0079(1) 


0.0080 


0.0080(1) 


5 


-0.0060 


-0.0059(4) 


-0.0060 


-0.0060(2) 



Having thus gained confidence in the validity of our IMC 
approach and its implementation we can move on to apply it 
to our objective namely to determine effective actions repro- 
ducing Yang-Mills thermodynamics, in particular the decon- 
finement phase transition. The first question we want to con- 
sider is whether the effective couplings Ai viewed as functions 
of the Wilson coupling (3 are sensitive to the phase transition. 
The answer turns out to be affirmative: Fig. 1 clearly shows 



0.15 




FIG. 1 : Behavior of the couplings Ai, . . . , As (appearing in the NN 
A-action) as a function of the Wilson coupling j3. (IMC based on 
algebraic SDEs.) 
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a rather drastic change in the behavior of the couplings at a 
value of j3\ ~ 5.69 for all Aj, i = 1, . . . , 5. Across the tran- 
sition, i.e. in both phases the dominant coupling is the 'funda- 
mental' one, Ai, followed by the octet couplings A3 and A5. 
The latter is actually a 'potential' coupling in the sense of [24] 
as it multiplies the center symmetric single-site octet charac- 
ter xn ( see App. A). The couplings A2 and A4 are clearly 
subdominant. 




0.0 1 1 : 1 1 1 

5.40 5.45 5.50 5.55 5.60 5.65 5.70 
P 

FIG. 2: Comparison between Yang-Mills and effective Polyakov 
loops for a lattice of size 8x3 and five NN couplings Ai, . . . , A5. 

The natural observable to address is, of course, the 
Polyakov loop which serves as the order parameter of the first- 
order SU (3) phase transition. In Fig. 2 we compare the ef- 
fective and Yang-Mills Polyakov loops for a relatively small 
lattice of size 8 3 x 3 where the would-be discontinuities (in in- 
finite volume) of the transition are still fairly smooth. Some- 
what surprisingly it is only the algebraic SDEs which repro- 
duce the behaviour of the Yang-Mills Polyakov loop reason- 
ably well. The geometrical SDEs, on the other hand, fail to do 
so, at least in the region just above the transition point. This 
is a first hint that there is some inherent instability in the IMC 
procedure - in particular if the geometrical SDEs are used. 

If we move on to larger lattices where the jump of the order 
parameter at the critical coupling becomes more pronounced 
one finds the behavior displayed in Fig. 3. Again, the alge- 
braic SDEs work satisfactorily unlike the geometric ones for 
which, in particular, the sudden rise of the order parameter 
appears at a larger value of (3, namely f3 gco ~ 6.5. This 
is substantially larger than the critical coupling, (3 C ~ 5.7. 
Our explanation for this behavior is the fact that, due to the 
first-order nature of the transition, there are rather sharp phase 
boundaries in the space of coupling constants, Ai. Hence, a 
tiny change in the couplings (presumably well within the IMC 
error bars) may easily have a large effect: by straying into 
the 'wrong' phase the Polyakov loop will suddenly explode 
or collapse. This nonlinear effect is rather difficult to evade 
considering the unavoidable (if small) instabilities of the IMC 
procedure. As a result, as we inevitably increase these inac- 
curacies by adding more coupling we expect this fine-tuning 



1.0 - 

0.8 - 
(W\) 




0.0 1 • 1 1 • 1 1 

4.5 5.0 5.5 6.0 6.5 7.0 7.5 8.0 

P 

FIG. 3: Comparison between Yang-Mills and effective Polyakov 
loops for a lattice of size 16 3 x 4 and five NN couplings Ai, . . . A5. 



TABLE II: IMC consistency check for the extended NN and NNN 
A-actions including up to 11 couplings A;. 





extended NN 


NN + NNN 


i 


Aj, input AiJMC 


Ai, input A;, IMC 



1 

2 
3 
4 
5 
6 
7 
8 
9 

10 
11 
12 
13 
14 
15 
16 



problem to become enhanced even further. The following sub- 
section will precisely address this topic. 

B. NN and NNN couplings 

There are (at least) two possibilities to generalise the NN 
A-action of the previous subsection. One may either extend 
the NN terms to higher group representations or include inter- 
actions of larger range, say NNN. 

The new NN terms we will add are of strong-coupling order 
P 3Nt , the additional NNN terms of order (3 2Nt . As we are 



-0.0050 


-0.0050(1) 


-0.0400 


-0.0400(1) 


0.0100 


0.0100(1) 


0.0100 


0.0099(1) 


-0.0150 


-0.0151(2) 


-0.0200 


-0.0201(1) 


0.0070 


0.0070(1) 


0.0300 


0.0300(1) 


-0.0080 


-0.0080(5) 


0.0050 


0.0052(3) 


0.0090 


0.0091(1) 






0.0030 


0.0030(1) 






-0.0030 


-0.0030(1) 






0.0080 


0.0081(1) 






-0.0060 


-0.0058(1) 






0.0020 


0.0020(2) 










-0.0020 


-0.0021(1) 






0.0060 


0.0060(1) 






-0.0030 


-0.0030(1) 






0.0040 


0.0039(1) 






-0.0070 


-0.0070(1) 
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working at large j3 one cannot predict which ones are going to 
be more important. Rather, this will be one of the questions to 
be considered in what follows. 

As before we have first tested the consistency of our SDEs. 
Table II shows once again that even for a total of the order of 
ten couplings the method works well: IMC output reproduces 
input for the effective theory. The empty input entries in Table 
II correspond to vanishing couplings. For a few sample cou- 
plings we have checked that vanishing input correctly entails 
vanishing output as well. 

To determine the behavior of the effective couplings Aj, 
i = 1, . . . , 16, as a function of the Wilson coupling (3 we have 
used a set of 4 x 10 6 configurations per (3 on a 16 3 x 4-lattice. 
This amounts to 5 x 10 4 (10 3 ) uncorrelated configurations far 
away from (close to) the phase transition. The IMC results 
are shown in Figs. 4 and 5 displaying the effective couplings 
as functions of (3. Again we note that the fundamental and 
octet potential couplings (Ai and A5) dominate in size. More 
important from a principal point of view, however, is the ob- 
servation that the behavior of coupling constants is very sen- 
sitive to the choice of operators. Let us compare, for instance, 
the coupling A3 in the two Figs. 4 and 5. For the extended NN 
ansatz (Fig. 4) it is comparable in magnitude with Ai and A5 
and behaves similarly as for the simple NN ansatz of Fig. 1. 
However, as soon as we include NNN operators its magnitude 
drops by 100% and its behavior changes drastically (Fig. 5). 
The latter is also true quite significantly for the octet coupling 
A5. This clearly signals an instability of the IMC methods, at 
least as far as the determination of the couplings is concerned. 



0.20 




-0.25 



FIG. 4: Behavior of the couplings Ai, . . . , An (appearing in the ex- 
tended NN A-action) as a function of the Wilson coupling /3. (IMC 
based on algebraic SDEs.) 

Nevertheless, it might still be possible that largely differ- 
ent sets of couplings lead to more or less identical behavior 
of observables. Comparing the behavior of the Polyakov loop 
in the effective and Yang-Mills theories rules out this possi- 
bility. As Fig. 6 shows the Polyakov loop when calculated 
in the effective models is extremely sensitive to the choice of 




FIG. 5: Behavior of the couplings Ai , . . . , A5, A12, • • • , Ai6 (appear- 
ing in the NN+NNN A-action) as a function of the Wilson coupling 
/?. (IMC based on algebraic SDEs.) 



operators and the value of (3 around (3 C . For both choices of 
SDEs the effective order parameter significantly overshoots 
the Yang-Mills one in a small /3-range near (3 C (see the spikes 
in Fig. 6). This means that in the space of effective couplings 
the phase boundary to the deconfined phase have slightly (and 
for a short range of (3 values) been crossed albeit with dras- 
tic effect due to the discontinuous behavior of the Polyakov 
loop. We conclude that the fine tuning problem encountered 
in the previous subsection indeed becomes more severe if we 
include more operators (and hence increase the instabilities of 
the IMC procedure). For the given number of effective cou- 
plings (of order ten) we have not been able to get this problem 
under control. 



3.0 r 



2.5 



2.0 



1.5 



1.0 ■ 



0.5 ■ 



0.Q L 



■ Yang-Mills 

-NNO(/3 M ') 
-NNN 




4.5 5.0 5.5 6.0 6.5 7.0 7.5 8.0 

P 

FIG. 6: Comparison between Yang-Mills and effective Polyakov 
loops for the extended NN-action (couplings Ai, . . . , An) and the 
NN+NNN action (couplings Ai, . . . , A5, A12, . . . , Ai 6 ). 
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V. SUMMARY AND OUTLOOK 

In this paper we have applied the IMC method to study the 
finite temperature phase transition of SU(3) Yang-Mills the- 
ory. A crucial input were novel Schwinger-Dyson equations 
based on algebraic and geometrical properties of the SU (3) 
Haar measure. The resulting equations constitute overdeter- 
mined linear systems for the effective couplings A, which 
were solved numerically via least square techniques. The 
method works well if the number of couplings is sufficiently 
small, say of the order of five. However, already in this 
case one notes a fine-tuning problem as the behavior of the 
Polyakov loop depends in a very sensitive and nonlinear way 
on the effective couplings. This fact can be traced back to 
the discontinuities associated with the first-order character of 
the phase transition. Hence, the problem becomes more pro- 
nounced in larger volumes. 

If we increase the number of effective couplings and thus, 
inevitably, the instabilities in their IMC determination, the 
fine-tuning problem again becomes more severe. This holds to 
such an extent that we could no longer gain numerical control 
and hence could no longer reproduce the Yang-Mills behav- 
ior of the Polyakov loop in the vicinity of the critical Wilson 
coupling, /3 = (3 C . We believe that an improvement on this 
situation will require nontrivial modifications of the IMC pro- 
cedure like e.g. smoothening of the loop in order to avoid the 
unphysical spikes of Fig. 6. In addition, it would be interest- 
ing to check whether Creutz's microcanonical demon method 
[15] mentioned in the introduction yields better results. 

We conclude, nevertheless, with the positive statement that 
the IMC method does work for the first-order £[7(3) transi- 
tion as well if we allow for only a small number of terms in 
the effective Polyakov loop actions. 



S*=J2 (xio(V x ) X 2o(Vy) 

+ X2o(V x )xio(V y ) + c.c.) , 

s 6 = $>ii (p«)> 



c.c. 



{xy) 

S-=J2 (X2l(P«)Xl2(P») + C.C.) , 

5 8 =E (X3o(P.)Xn(P«) 

{xv) +Xu(P«)X3o(P»)+c.c), 

S 9 = (X2l(V x )X20(Vy) 

{XV) + X20(V x ) X 2l(Vy) + c.c.) , 

Sio=J2 (X2l(Px)X0i(Py) 

{XV) + Xoi(V x ) X 2i(Vy) + c.c.) , 
Sn=Y,(X3o(V x )+c.c.) , 

X 

Si2 = Y,(xio(r x )xoi(r z ) + c.c.) , 

[SBZ] 

Sis = (xw(V x )xoi(V z ) + c.c.)xu(Vy) , 

{ X yz) 

Su= (xw{V x )xm{V v )xw{V z ) + c.c.) , 
{xyz) 

Sis = Y (Xio(V x )xw(V y )xio(V z ) + c.c.) , 

{xyz) 

Sie= J2 (xio(V x )xoi(V y ) + c.c) 

{xv ' vw) ■(xio(P«koi(^)+c.c.) 



(A4) 

(A5) 
(A6) 
(A7) 
(A8) 

(A9) 

(A10) 

(All) 
(A12) 

(A13) 

(A14) 

(A15) 

(A16) 
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APPENDIX A: OPERATORS FOR THE A-ACTIONS 

In this paper we have used up to 16 different operators ap- 
pearing in the A-action (7): 

51 = Y.(xio(V x )xoi(V y )+c.c.) , (Al) 

(xy) 

5 2 = J2(X3o(V a )x(a(V v ) +c.c.) , (A2) 

(xy) 

S3=£>ii(7>.)Xii(Pv). (A3) 
(xy) 



The NN and NNN relationships are denoted in terms of brack- 
ets the meaning of which is explained in Fig. 7. Hence, the 
operators Si to Su obviously describe (extended) NN inter- 
actions, while S12 to Si g are NNN terms. In a strong coupling 
(small-/?) expansion the terms Si, ■ . ■ , S5 and S12, . . . , Sw 
would be of C(/3 2Ar «), the terms S 6 ,...,S U of <D(f3 3Nt ) [31]. 



APPENDIX B: OPERATORS FOR THE k- ACTIONS 

If we extend the strong coupling NN contributions to 
0((3 3Nt ) the effective action becomes a series of nine terms, 

9 

S e s — n a I a , (Bl) 

a=l 

which is referred to as the extended NN action (as is its A- 
equivalent, see Appendix C below). 

If we allow for NNN interactions (which, however, do not 
extend beyond single plaquettes) up to order 0((3 2Nt ) we end 
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1111 



x w 



y v 
* • <> <> 



I I I I I 



FIG. 7: The neighboring relationships of the marked sites corre- 
spond to the bracket notations [xv], [yw], (xyv), (yvw), (vwx), 
(wxy), (xy, vw) and (xw, yv). 



ho= (Xio(P-)xoi(P«) + c.c) (B12) 
{XVZ) -(xio^koiC^+cc.), 

Ai= I] (xio(Px)xoi(^) + cc.) (B13) 
{XV ' VW) •(xio(K)xoi(^)+c.c.). 



APPENDIX C: LINEAR COUPLING RELATIONS 



For the extended NN action (Bl) the relation between the 
\ and the n a is 



up with the effective action 



Sett — ^ ' K ala 

oe{l,2,3,4,10,ll} 



(B2) 



This (and its A-equivalent, see Appendix C below) is referred 
to as the NN+NNN action. 

The resulting operators are given by 

h = J2 (Xio{V x )xoi{V v ) + c.c.) , (B3) 

(xy) 

h = J2 (X2o(^)xo 2 (^ y ) + cc.) , (B4) 
(xy) 

h = Y,Xu(V a ,)xii(V v ), (B5) 
(^y) 

h = J2(xio(Px)xoi(Py) +cc.) 2 , (B6) 

(xy) 

h = Y< (X3o(V x )xo3{V y ) + c.c.) , (B7) 

(xy) 

h = Y.(X2i(V x )xi2(V y ) +c.c.) , (B8) 

(xy) 

h = Yl (Xio(r x )xoi(Vy) + c.c.) (B9) 

{XV) ■ (X20(V x ) X 02(Vy) + c.c.) , 

h = (Xio(^)xoi(^) + c.c.) (BIO) 

{XV) ■Xll{V x )Xll{Vy), 

h = ^(xio(^)xoi(^)+cc.) 3 , (Bll) 
(xy) 



(Xi) 




(l 
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The analogous relation for the NN+NNN action (B2) is 
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